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Multilayer fractional quantum Hall wave functions can be used to construct the non-Abelian 
states of the Z& Read-Rezayi series upon symmetrization over the layer index. Unfortunately, this 
construction does not yield the complete set of Z& ground states on the torus. We develop an alter¬ 
native projective construction of Z& Read-Rezayi states that complements the existing one. On the 
multi-layer torus geometry, our construction consists of introducing twisted boundary conditions 
connecting the layers before performing the symmetrization. We give a comprehensive account of 
this construction for bosonic states, and numerically show that the full ground state and quasihole 
manifolds are recovered for all computationally accessible system sizes. Furthermore, we analyze 
the neutral excitation modes above the Moore-Read on the torus through an extensive exact di- 
agonalization study. We show numerically that our construction can be used to obtain excellent 
approximations to these modes. Finally, we extend the new symmetrization scheme to the plane 
and sphere geometries. 


I. INTRODUCTION 

Exotic correlated and topologically ordered quantum 
states can often be constructed from less exotic parent 
states via the action of specific many-body projection 
operators. Projective constructions can yield topological 
order by starting with entirely uncorrelated single Slater 
determinant parent states of (topological) band insula¬ 
tors. Well-studied examples are so-called parton com 
structions of fractional quantum Hall (FQH) statef® 
and more generally fractional topological insulator^ 7 9 . 
Other projective constructions can change the form of 
topological order in a parent state to a more exotic one, 
for example by turning an Abelian FQH state into a 
non-Abelian one. This has been exploited in topologi¬ 
cal orders described by conformal field theories via the 
so-called coset projections 10 13 . 

The Zfc Read-Rezayi seried^ is a well known sequence 
of FQH states with non-Abelian topological order. The 
k = 1 member of this series and parent state for our con¬ 
struction is the (Abelian) Laughlin statd^. The k = 2 
member of this series is the Moore-ReacP^ state support¬ 
ing Majorana excitations and a well studied candidate 
state for the FQH plateau of electrons at v = 5/2. The 
k = 3 member of this series supports Fibonacci anyons, 
which can in principle be used to perform the opera tions 
of a universal topological quantum computer 1718 . For 
bosons, the Z/~ Read-Rezayi series can be obtained using 
projective constructions, by starting from the Laughlin 
state as a parent state. Reference [12] considered k in¬ 
dependent layers of Laughlin states. The projective con¬ 
struction then consists of symmetrizing over the layer de¬ 
gree of freedom, yielding a single layered Z& Read-Rezayi 
state. The approach of Ref. 12 was based on conformal 
field theory arguments, and was subsequently tested by 


numerical simulations on the sphere geometry 1 ^. Besides 
the exact zero-energy ground states and quasihole states, 
this procedure also provides accurate trial wave functions 
for the low energy neutral and quasielectron excitations 
of the Zfc Read-Rezayi states on the sphere geometry, 
as numerically shown in Refs. I20ll23l More than a use¬ 
ful mathematical trick, the projective constructions may 
also be of physical importance: several proposals®^ to 
build non-Abelian order from Abelian systems are indeed 
based on this construction, using either a tunneling term 
or an interaction between layers. 

However, this established multi-layer symmetrization 
procedure is not capable of constructing the Z& Read- 
Rezayi states for all system sizes on the torus geometry. 
For instance, the Moore-Read ground state appearing for 
an odd number of particles on the torus cannot be writ¬ 
ten using this procedure, and similarly, other states in the 
Zfc series fail to have a multilayer symmetrization descrip¬ 
tion. In this paper, we present an alternative projective 
construction that allows us to also obtain those miss¬ 
ing ground states. We consider a /^-layered torus with 
twisted instead of periodic boundary conditions, which 
can be realized through an extended topological defect 
that connects the layers. In addition, the construction 
can be used to obtain the quasihole states of the Z& Read- 
Rezayi series for all system sizes on the torus. We show 
that the multilayer torus with twisted boundary condi¬ 
tions is equivalent to an enlarged, single-layer torus. Our 
construction thus enables us to extract any bosonic Read- 
Rezayi state from a single Laughlin wave function. We 
generalize this novel scheme to the plane and sphere ge¬ 
ometries as an alternative to the symmetrization over k 
independent layers. 

A second part of our study is dedicated to the neutral 
excitation modes above Z& Read-Rezayi states. These 
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mode s are the higher-/^ cousins of the magnetoroton exci¬ 
tation^ 4 23 above the Zi Laughlin state. We detail their 
structure on the torus geometry for the simplest non- 
Abelian case, the Moore-Read state. Furthermore, we 
show that symmetrization of the magnetoroton excita¬ 
tions above the Laughlin state yields excellent approxi¬ 
mations to the excitation modes above the Moore-Read 
state both on the torus and the sphere. While these trial 
wave functions are different from those obtained using 
the multilayer approactPHHSS they describe the neutral 
mode with the same accuracy. 

This paper is structured into three parts. Sect ion |ll|de- 
fines the projective construction on the torus with twisted 
boundary conditions, interprets it both geometrically and 
through its action on the orbitals, and presents the ex¬ 
act numerical results on zero energy states. The subse¬ 
quent Sec. [HI] is concerned with the neutral excitation 
modes above the Moore-Read state on the torus, that 
are accessed both with exact diagonalization and with 
the projective construction. Finally, Sec. |IV| contains a 
detailed discussion of the results on the plane and sphere 
geometries, including both the zero-energy states and the 
neutral excitation modes. 


II. READ-REZAYI STATES ON THE TORUS 
A. Ground states and their degeneracy 

The bosonic Z& Read-Rezayi® state |T/~) is the dens¬ 
est zero-energy ground state of the (fc + l)-body contact 
interaction 

flint,k = J d 2 r : [p(r)] k+1 (1) 

when projected into the lowest Landau level. Here, 

p(r) = ip\r)ip(r), (2) 

is the density operator written in terms of the field op¬ 
erator r ), that creates a boson at position r, and :: 
represents normal ordering. 

On the torus, the Z& Read-Rezayi ground state appears 
at filling v = k/ 2, that is, if the number of particles N 
and the number of flux quanta obey 

kN. 0 = 2 N. (3) 

This relation cannot be satisfied for all system sizes, since 
both A and have to be integer. For example, if k is 
odd and N ^ is odd, there exists no ground state. If a 
ground state exists, it can exhibit a topological degener¬ 
acy on the torus (but not on the sphere). This degeneracy 
is a hallmark of the specific topological order and is in 
one-to-one correspondence with the number of topologi¬ 
cal anyon excitations of the phase. The topological order 
realized by a bosonic Z& Read-Rezayi state is labelled by 
the affine Lie algebra SU(2)/ C and contains (k + 1) irre¬ 
ducible representations, corresponding to (fc + 1) topolog¬ 
ical anyon excitations. Thus, on the torus with an even 



k odd 

(ksjkj,) deg. 

k even 

(ks,k y ) deg. 

N(j> odd 

- 

(0,0) 1 

even 

(0,0) *±i 

(0,NJ2) *±i 

(0,0) f 

(0,AW2) 1 

(AV/2,0) 1 

(JV2.JV2) |-1 


TABLE I. Ground state degeneracies and momenta of the 
Zfc Read-Rezayi state on the torus. For k odd (respectively 
even), the momentum quantum numbers k x , k y thus belong to 
a N#/ 2 x N<f> (respectively N$ x N#) Brillouin zone. These mo¬ 
mentum quantum numbers are defined in the Landau a>gauge 
with vector potential A(r) = (0 ,—Bx). See Appendix [A] for 
the definition of the momentum quantum numbers. 


number of flux quanta 7V^, any Z& Read-Rezayi state is 
(k + l)-fold degenerate. The momentum quantum num¬ 
bers of the degenerate ground states are given in Tab. [I] 
and Appendix [A] Note that there is a fundamental dif¬ 
ference between k even and k odd. For k odd, N must 
be a multiple of k for the filling fraction to be v = kj 2. 
When k is even, there are two alternatives: N is either 
an integer multiple of k (and N<f, is even) or a half-integer 
multiple of k (and is odd). In the latter case, there 
is a unique Z& state. Note that this feature is specific to 
the torus geometry, and absent on the sphere. 

If one deviates from the filling factor v — k/2 by in¬ 
creasing the number of flux quanta, quasihole excitations 
are nucleated. Their wave functions correspond to the 
zero energy states of the model Hamiltonian Q. The 
number of zero energy states (and their quantum num¬ 
bers) for a given number of particles and a given number 
of flux quanta can be deduced from the clustering prop¬ 
erties of these stated 6 37 1 (for the sphere geometry, closed 
formulas are known 38 40 ). 


B. Zfc Read-Rezayi states from symmetrization 

In this work, we exploit the fact that a 7L g Read- 
Rezayi state I'kgk) can be obtained by symmetrizing over 
g independent copies of a Z& Read-Rezayi state IT&). 
This procedure was first introduced for the Moore-Read 
state, which can be written as a bilayer Laughlin state 
in the strongly paired regimd- 41 42 1 In the conformal 
field theory language, the coset construction described 
in Refs. uni and [IT] (and later generalized to the Read- 
Rezayi series in Ref. mu uses Abelian theories to build 
the parafermion Hall states. 

Introducing a symmetrization operator S g ^ i, we write 

i« 9fe ) = 5 9 _ >1 (i«fe)®---®i$ fc >y (4) 

1 V ^ 

g times 

Here, we can picture the direct product of |T&) states as 
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being defined on g “layers” of the same manifold. For 
example, we can construct a 7L g Read-Rezayi state for 
any g by symmetrizing over g layers of v = 1/2 bosonic 
Laughlin states. 

Let be any basis of single-particle bosonic creation 
operators acting on the single particle Hilbert space Ht 
in the position basis on the torus and denote by (j^- v 
l = 0, • • • ,g — 1, the corresponding operators belonging 
the Z-th layer of Ht in the direct sum Ht ® • • • ® Ht 
that constitutes the single-particle Hilbert space of the g- 
layered torus. Then, the symmetrization S g ^i is a map 
from Ht ® • ■ ■ ® Ht to Ht defined by 

4>),i Sg ^ 1 4>\- (5) 

We note that there is a fundamental problem with the 
symmetrization construction Q for systems of finite size. 
If a densest ground state l^/Jhas a particle numbers N 
that is not divisible by g , this particle number cannot 
be spread equally over g densest states | T&). To give a 
simple example, consider the bosonic Moore-Read state 
|^2) with N = 7 particles, which is a zero-energy ground 
state of the Hamiltonian 0 on a torus with = 7 flux 
quanta (filling v = 1). On the same torus, the bosonic 
v = 1/2 Laughlin state |Ti) does not exist as a zero en¬ 
ergy ground state, since = 7 is not divisible by two. 
One could consider symmetrizing over a Laughlin quasi¬ 
hole state with N = 3 and a Laughlin quasielectron state 
with N = 4 (the latter being at finite energy). However, 
this construction can at best yield an approximation of 
the desired ^2) ground state (we checked numerically 
that the corresponding overlap was of the order of 0.998 
for TV 1= 17). It might be possible to obtain the exact 
state by using the composite fermion^ expression of the 
quasiele ctron states, but these are very hard to obtain on 
the torn J 44 * 4 ^ . 

A similar obstruction appears for all Z/~ Read-Rezayi 
states with k even on a torus with odd N^. Adding flux 
quanta to the original v = k/2 system is a way to nucle¬ 
ate quasihole excitations of the Z/~ Read-Rezayi states. 
While the obstruction disappears for a sufficient number 
of added flux quanta, not all quasihole states are im¬ 
mune to this issue. The system with 7 particles on a 
torus pierced by = 5 flux quanta has a filling frac¬ 
tion strictly lower than 3/2, and thus admits zero energy 
states of the 4-body interaction of the Z3 Read-Rezayi 
state. However, the 7 particles can at best be spread 
into two layers with 2 particles (Laughlin quasihole state) 
and one layer with 3 particles (Laughlin quasielectron), 
resulting in a non-zero energy state after symmetrization. 
In general, we see that for odd and the largest N such 
that N/Nfi < k/2, the Z& states cannot be constructed 
using the symmetrization technique. 

A main result of this paper is an alternative sym¬ 
metrization scheme on the torus that remedies this ob¬ 
struction. We show that this scheme can also be used 
on the sphere as an equally powerful alternative to the 
multi-layer symmetrization. On the torus, the key idea 



FIG. 1. Bilayer system with different types of topological de¬ 
fects. a) Bilayer system with no defect, b) Bilayer system 
with a twist defect along the y direction. If the system has 
periodic boundary conditions, it is transformed into a mono- 
layer torus with a doubled length in the x direction, c) Bilayer 
system with a defect along the x direction. If the system has 
periodic boundary conditions, it is transformed into a mono- 
layer torus with a doubled length in the y direction. 


is to change the boundary conditions between the layers 
from periodic to twisted, as depicted in Fig. [l] Equiva¬ 
lently, this can be seen as introducing a topological de¬ 
fect line that permutes the layer indices^®*. Similarly 
to the symmetrization over multiple independent layers, 
our construction with twisted boundary conditions in one 
direction between the layers alone does not yield the com¬ 
plete Zfc manifold. To obtain the complete manifold, two 
of the three symmetrization schemes (untwisted multi¬ 
layer, twist in x direction, twist in y direction) have to 
be combined. 

As a manifold, a double layer torus of size L x x L y with 
twisted boundary conditions in the x-direction (respec¬ 
tively ^/-direction) is equivalent to a single layer torus 
of size 2 L x x L y (respectively L x x 2 L y ). The sym¬ 
metrization is then taken over the particle coordinates 
with 0 < x < L x and L x < x < 2 L x (respectively with 
0 < y < L y and L y < y < 2 L y ). This carries over 
to a g layered torus, for boundary conditions that fully 
permute the layer coordinates. (See Fig. 2] for an illus¬ 
tration of different possible boundary conditions of three 
layers.) In the next section, we will show that this equiv¬ 
alence between manifolds is also respected by the model 
Hamiltonian 0 for the Zfc Read-Rezayi states. More¬ 
over, we will work out the action of the symmetrization 
operator S g ^\ on the basis states of the g times larger 
tori: gL x x L y (the ( gT) x toms) and L x x gL y (the (gT) y 
torus). 

One may wonder whether yet another relevant sym¬ 
metrization scheme is obtained by imposing twisted 
boundary conditions in both x and y directions in a dou¬ 
ble layer system. As we illustrate in Appendix [Bj the re¬ 
sulting surface is equivalent to a single-layer toms whose 
spanning vectors form a given relative angle and norm 
ratio that depend on the symmetrization scheme that 
will be applied. The symmetrization schemes are then 
the same as for twisted boundary conditions in one di¬ 
rection, with the difference that the symmetrized state is 
also defined on a toms with that different angle between 
its spanning vectors. 
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(a) (b) (c) 



FIG. 2. Trilayer system with different types of topological 
defects, a) Trilayer system with no defect, b) Trilayer system 
with a defect along the y direction connecting the two up¬ 
per layers. With periodic boundary conditions, this system 
is equivalent to a bilayer torus, with the first layer having as¬ 
pect ratio 2, and the second layer having aspect ratio 1. c) 
Trilayer system with a defect along the y direction connect¬ 
ing respectively layers 1 (blue) and 2 (red), 2 and 3 (green), 
3 and 1. With periodic boundary conditions, this system is 
equivalent to a torus with aspect ratio 3. 


C. Equivalence between twisted boundary 
conditions and enlarged torus 

We want to show that the model Hamiltonian 0 
for Z/c Read-Rezayi states on the ^-layered torus with 
twisted boundary conditions is equivalent to the same 
Hamiltonian on a g times larger torus. For this, we mo¬ 
mentarily neglect the projection on the lowest Landau 
level, which will later be used to obtain all numerical 
results. We will see that the equivalence holds for the 
entire system and thus also for the lowest Landau level. 

Let us for concreteness fix the gauge to be the Lan¬ 
dau x-gauge with vector potential A(r ) = (0, —xB) and 
consider separately the cases of twisted boundary condi¬ 
tions in the ^/-direction and twisted boundary conditions 
in the x-direction. A complete set of single-particle op¬ 
erators on the g-layer torus T = [0, L x ) x [0,L y ) with 
(fully) twisted boundary conditions is given by the op¬ 
erators i/j] (r), that create a boson at position r G T in 
layer / = 0, ••• ,g- 1 . 

If twisted boundary conditions are applied along the 
x-direction, they obey 

( r + L x e x ) = e 2 ^ y/i ^ + 1 )mods (r), ( 6 a) 

where = BL x L y . A complete set of single-particle 
operators on the g times larger toms ( gT) x = [0 ,gL x ) x 
[0, L y ) is given by the operators ?/h(r), that create a bo¬ 
son at position r G ( gT) x . They obey periodic boundary 
conditions 


ip^(r + gL x e x ) = e 27rigN(t>y ^ Ly ( 6 b) 

If twisted boundary conditions are applied along the 
^/-direction, they obey 

i>\(r + L v e y ) = A+i)mod g (?)- ( 7a ) 

A complete set of single-particle operators on the g times 
larger torus (gT) y = [0 ,L X ) x [0 ,gL y ) is given by the op¬ 
erators V^( r ), that create a boson at position r G ( gT) y . 
They obey periodic boundary conditions 

ip\r + gL y e y )=ip\r). (7b) 


The model Hamiltonian for the Z& Read-Rezayi state 
on either manifold decomposes into a noninteracting Lan¬ 
dau level Hamiltonian and an interaction part 

Hk — Ho T H\ nt,fc) H}~ = Hq T -Hint,/c5 (8) 

that act on the Fock space build from the single-particle 
operators in Eq. (| 6 b|)/Eq. (pfb| ) and Eq. (| 6 a|)/Eq. ( |7a| ), re¬ 
spectively. Here, Hq and Hq are the single-particle opera¬ 
tors for the T and ( gT)i , i = x,y, tori, respectively. The 
respective interacting parts iLmt,/c and iAnt,/c are given 
by Eq. 0 and 



with the latter expressed in terms of the density operator 
pi(r) = in the Z-th layer. 

If the boundary conditions are twisted in the x- 
direction, the identification 

S fp\r + lL x e x ), 1=0, reT, 

( 10 a) 

provides a mapping between the two single-particle 
Hilbert spaces (and thus also the Fock spaces) under 
which iJint,*; is exactly mapped into H[ nt ^. Analogously, 
if the boundary conditions are twisted in the ^/-direction, 
the identification is 

^\{r) =^\r+ lL y e y ), 1 = 0, •••,0-1, reT. 

(iOb) 

The ultra-locality of the contact two-body interaction is 
crucial to make this connection to a g times as long torus 
for the interacting system. Under this same mapping, 
the single-particle Hamiltonian 




2 _ 


eAi(r) 


i(r ) 


with 

V ? A Ai(r) = B, 1= 0, • • • ,g - 1, reT 


maps into 

Hq = f d 2 r [iV r — eA(r)} 2 ^(r), 

J(gT)i 

with i = x, y, and 

A l (F) = A(r + lL l e i ) 


(ii) 


( 12 ) 


(13) 


(14) 


defining the gauge potential of (gT)i that thus also obeys 
V r A A(r) = B for every r G ( gT\. 

With this, we have shown the equivalence of the inter¬ 
acting model Read-Rezayi systems on the ^-layered torus 
with twisted boundary conditions and on the g times 
larger toms. Since the magnetic field B is preserved un¬ 
der this mapping, the g times larger toms is pierced by g 
times as many flux quanta, i.e., gN as the initial toms. 
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D. Symmetrization on the torus with twisted 
boundary conditions 


To implement the symmetrization operation, we want 
to compute the action of the operator on a 

many-body wave function defined on the g times larger 
torus ( gT)i , i = x,y [( gT\ is the torus enlarged g times 
in the i direction]. The action of can be written 

in first quantized notation. Let us denote the translation 
operator 72* e4 that translates single-particle operators 
by Liei on ( gT)i , i = x,y. The symmetrization identifies 
positions that are related by any magnetic translation 
TsLid , where 0 < S < g, S G Z. Starting from a many- 
body wave function T on a (gT)i torus, the symmetrized 
wave function reads 

(fi> • • • ,vn) = 

( N 

^ > 11 

{0<5 3 '< fl } V=1 



where rj denotes the coordinates of the j-th particle re¬ 
stricted to the L x x L y region and the j label on the 
translation operator indicates it acts on the j -th parti¬ 
cle. With our gauge choice, Eq. (15) has an even more 
explicit expression for S^ g T) y ^T 


with 


<t>\r) =e 27rijv ^/ i «^ t (f ■ 
0 f (r) =4>\r + L v e y ). 




(17b) 


Thus, these single particle operators can be defined on a 
g times smaller torus. Using the identification (10), we 
obtain the identity 


S(gT)i^T [^( r )] — S(gT)i^ T ( r + IL^i)] , 

v Z = 0 , - - - ,g- 1 , 


(18) 


for the operators V^( r ) on ( gT)i , i = x,y. Let us in¬ 
terpret this result in terms of the eigenspaces of the 
translation operator 72* e 4 - It has eigenvalues e 2lTls 
s = 0, • • • ,g — 1. We observe that S( gT ).^ T is precisely 
the projector on the single-particle eigenstates of 72^ 
with eigenvalue 1 on ( gT)i . (Note that the projection 
on any eigenvalue e 2?rls is an equally good option for 
the symmetrization as we will discuss below). While this 
result has been derived for the entire Hilbert space, it 
is also true in any 72 i er^ nvar ^ an t subspace, such as the 
lowest Landau level. 

In the x-Landau gauge, a basis </>J(r) of the lowest 
Landau level on the torus T with flux quanta is given 
by 


S ( gT) y ^ T y (r ir -- ,r N ) = 

Y ’J' (n + 5iL y e y ,- ■ ■ ,r N + 5 N L y e y ) ( 16 ) 

{0<5 i<s } 

The vanishing properties of the symmetrized wave func¬ 
tion can immediately be deduced from those of T. In¬ 
deed, if T is a Laughlin wave function (or one of its quasi¬ 
hole excitations), it vanishes when two particles are at 
the same point. There are terms in the sum of Eq. (16) 
where , S g are all distinct. If we put the g first 

particles at the same position, these terms do not need 
to vanish. We now bring the (g + l)-th first particles 
at the same position. In each term of the sum, at least 
two Si among 5i, • * • , 5 g +1 are equal, ensuring that T is 
evaluated at a position including two equal coordinates. 
Thus S( g T) y ^T^ vanishes when g + 1 particles are at the 
same position, proving it is a zero energy eigenstate of 
the (g + l)-body contact interaction ([!]), i.e. a 7L g Read- 
Rezayi wave function (or one of its quasihole excitations). 
A similar argument holds true for the symmetrization on 
the ( gT) x torus, including Eq. ( p~5] ) . 

A more practical expression of these symmetrization 
operators can be obtained by considering only their ac¬ 
tion on the single-particle operators or basis states. Using 
the definition of S g ^i from Eq. and the boundary 
conditions (6a) and f7a|), we see that indeed 


maps into operators $(r) that obey the correct bound¬ 
ary conditions on T, irrespective of whether the initial 
torus is (gT) x or ( gT) y , namely 








VZ = (),•■■ ,5-1, ( 17a ) 


^( p )~*7FT ^ 

V L yVir m =-oo 


e iy(27T j/Ly+mLx) 


x e 


— (cc+27T j/Ly+mL x y/2 


(19) 


with the orbital quantum number j = 0 , • • • , — 1 . 

(We will consider single-particle states instead of second 
quantized operators from here on.) The basis (j)f T ^ x (r) 
of the lowest Landau level on (gT) x with gN^ flux quanta 
is obtained from Eq. (19) by replacing L x with gL x and 


thus j may take values j = 0, • • • , gN^ — 1. The action 
of Tl x & x on the basis states of ( gT) x is given by 


TL X 


k (gT)a 


l(9 t )x 

'(j+N^modgNtf, 


( 20 ) 


A basis for the eigenspace of Tl x & x of eigenvalue 1 in 
the lowest Landau level is given by the equal-amplitude 
superposition of basis states with momenta j that are 
apart. With the help of Eq. (19) one can explicitly check 
that 



1 


9~ 1 



j = 0, • • • , — 1. (21) 


The interpretation of this result is that symmetrization 
of a many-body state over the orbitals that are spaced 
on the ( gT) x torus is equivalent to symmetrizing over 
the layers of the ^-layered T torus with twisted boundary 
conditions in the x direction. 

Next, we want to study the case where the twist in 
the boundary conditions is along the y direction, while 














6 



FIG. 3. Schematic representation of the symmetrization pro¬ 
cedure to obtain the ground state on the T torus of dimensions 
L x x L y pierced by N# flux quanta. The orbitals colored in 
red have become equivalent after symmetrization. a) The T 
torus bilayer of lengths L x x L y , pierced by N, 0 flux quanta, 
b) The (2 T) x toms of lengths 2 L x x L y , pierced by 2flux 
quanta, c) The (2 T) y torus of lengths L x x 2 L y , pierced by 
2 N 0 flux quanta. 


keeping the Landau gauge fixed in the x-direction. We 
are now seeking the eigenspace of TL y e y on the (gT) y 
torus. From Eq. (19) we find 


7L 


dgT) t 


0 2nij/g 




( 22 ) 


include this optional parameter 
S(gT)i^T^ ( r h • ’ • 5 r N ) = 

(J^L 2ni8 as \ 

II e ,r N )- 

{0<5j <g} \j= 1 / 

(26) 

We note that setting s is akin to projecting on the 

27ris 

eigenspace of Th iei with Tl iSi eigenvalue e » on (gT)i. 


E. Numerical results on zero-energy states 

We will now discuss the exact numerical results on 
the symmetrization construction in detail for the ground 
states and quasihole states of the Read-Rezayi series on 
the torus, i.e., for all zero energy states of the model inter¬ 
action 0- While we have shown in the previous section 
that the symmetrization leads to Read-Rezayi states, we 
still have to investigate in which cases the procedure al¬ 
lows to recover all the states. 


Thus, the eigenvalue 1 subspace of TL y e y is spanned by 
the basis states with j = jg , where j = 0, • • • , — 1. In 

other words, the states 




j = 0 , • • • , Afy - 1 


(23) 


form a basis for the T torus. The interpretation of this 
result is that the projection of a many-body state on the 
orbitals that are spaced g on the (gT) y torus is equivalent 
to symmetrizing over the layers of the ^-layered T torus 
with twisted boundary conditions in the y direction. 

Returning to second quantized formulation, we can 
summarize the action of S( 9 t) x ^t as follows. On the 
torus ( : gT) x = [0 ,gL x ) x [0, L y ) 




"j mod 5 


(24) 


i.e., a many-body state has to be symmetrized over all 
orbitals with quantum numbers j differing by N^. On 
the torus (gT) y = [0 ,L X ) x [0 ,gL y ) 




j-s)/g ’ 

0 , 


j mod g = s, 
else, 


(25) 


i.e., a many-body state is projected to the orbitals with 
quantum numbers j that are multiples of g (for 8 = 0 ). 
In choosing different 8 = 0,••• , <7 — 1, we obtain a family 
of equivalent symmetrization operators which only differ 
by a global translation between the origin of the coor¬ 
dinate systems on T and ( gT ) y . These symmetrization 
schemes are schematically illustrated in Fig. [3] Note that 
this freedom in the choice of s has an analogue in first 
quantized notation. Indeed, we can rewrite Eq. (15) to 


1. Symmetrization with periodic boundary conditions 


We start by testing the previously known 12 multi-layer 
symmetrization construction for the case of two layers of 
Laughlin ground states on the torus. We choose a system 
with an even number of particles, so that they can be 
evenly distributed among the two layers. 




( k y + k y) mod N 4> 
2 


^ — <?2^1 




(27) 


Notice that the states carry center of mass momentum 
quantum numbers k y and k ' y , which can take the values 0 
and N^/2 for the two degenerate Laughlin states. As the 
symmetrization makes the layers indistinguishable, the 
choices (k^, k^) = (0,A^/2) and (k y ,k! y ) = (A^/2,0) de¬ 
liver the identical Moore-Read state. In contrast, the 
choices {k y ,k' y ) = (0,0) and (k y ,k' y ) = (TV^/2, A^/2) 
yield two different Moore-Read states in the same k y 
sector. The latter two states are not eigenstates of the 
relative translation operator 7^ rel , even if the Laughlin 
states were (see AppendixjA]) . The diagonalization of 7^ rel 
in the subspace defined by the two degenerate Moore- 
Read states yields one state with k x = 0 and one with 
k x = N ( f ) /2. 

Going beyond the ground states, we numerically 
checked that one can obtain a Moore-Read quasihole 
state by symmetrizing two decoupled systems of Laugh¬ 
lin 1/2 quasiholes (i.e. one flux quantum added com¬ 
pared to the ground state filling fraction). As this pro¬ 
cess can sometimes be redundant, we have to extract a 
linearly independent basis from all states obtained after 
the symmetrization. We compare the number of linearly 
independent symmetrized products of Laughlin quasi¬ 
holes states with the number of Moore-Read quasihole 
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states for different values of TV < N^. For each < 12, 
we considered the systems with 4 < N < N^. For all 
even N < TV^, and all odd N < 7V^, we found that 
all Moore-Read quasihole states can be constructed as a 
symmetrized product of Laughlin quasihole states. 

As a result, the only zero-energy state of the three- 
body contact interaction that cannot be reproduced us¬ 
ing the symmetrization construction over two indepen¬ 
dent layers is the Moore-Read ground state that lies in 


the N = odd sector that we discussed in Sec. IIB 


2. Symmetrization with twisted boundary conditions 


We will now numerically test the symmetrization con¬ 
structions on the g-layer torus with twisted boundary 
con dition s, or equivalently the g times as long torus (see 
Sec. IID). We will start with a Z& state on the ^-layered 


torus to construct a Z g k state on the single-layer torus 
by symmetrization. Not only will this construction al¬ 
low us to access the missing Z g k states for odd, it 
also yields the entire quasihole manifold of Z& states. 
This construction comes with a trade-off: some of the 
Z g k ground states accessible using the multilayer (with 
periodic boundary conditions) representation cannot be 
reached using only one type of twisted boundary condi¬ 
tions. It is simply the consequence of the symmetrization 
operator yielding one state per vector it acts on. For in¬ 
stance, the Moore-Read state is constructed by choosing 
g = 2, k = 1. Acting with S^l x ^l x on the twofold degen¬ 
erate Laughlin state defined on the ( gT) x torus will yield 
both Moore-Read states in the k y = 0 sector, but not the 
Moore-Read state in the k y = N/2 sector. Conversely, 
only the Moore-Read states in the k x = 0 sector can 
be constructed by acting with S^Ly^Ly on the Laughlin 
state defined on the (gT) y torus. In summary, all of the 
Z g k states can be constructed using the symmetrization 
method, as long as at least two out of the three bound¬ 
ary conditions (no twist, twist in x direction, twist in y 
direction) are used. 

Enlarging the torus length in the x-direction — The 
numerical implementaion of S^ q t) x ^t is described in Ap¬ 
pendix^ This symmetrization scheme conserves the cen¬ 
ter of mass momentum k y modulo A^, but does not con¬ 
serve the relative momentum k x of a many body state 
(see Appendix [P] for details). 

We checked numerically that the application of 
S( gT } x ^T to the degenerate Z& state defined on the (gT) x 
torus yields only zero energy states of the (gk + l)-body 
model interaction. Focusing first on the densest Z& state, 
we apply S( 9 t) x ^t to the Laughlin state. We numerically 
checked the above property for N < 14 (g = 2), N < 12 
(g = 3), N < 14 (g = 4). This implies that applying 
<$2 l x ^l x to the Moore-Read state defined on the (2 T) x 
torus yields part (respectively all) of the Z 4 ground state 
manifold when is even (respectively odd). 

Decreasing the filling fraction, we have checked that 
the above property also holds true for quasihole states. 


In this latter case, we observe an additional property: 
the completeness of the symmetrized quasihole manifold. 
In other words, any zero energy state of the (gk + 1)- 
body interaction with N and such that N/N^ < gk/ 2 
can be reached by applying S^ 9 t) x ^t to the subspace of 
zero energy states of the (k + l)-body interaction, with 
the same number of particles and gN^ flux quanta. We 
numerically verified this statement in a number of cases 
for k = 1 and g = 2, 3, 4, including all cases with N <7 
and up to 3 added flux quanta (compared to the number 
of flux quanta for the densest Z/^ state). For gk = 3, 4, 
we also numerically verified this property for N < 10 and 
the smallest fraction of flux added to obtain a quasihole 
state. 

Enlarging the torus length in the y- direction — The 
numerical implementaion of S^ g T) y ^T is described in Ap¬ 
pendix [C] This symmetrization scheme does not conserve 
the center of mass momentum k y , but conserves the rela¬ 
tive momentum k x of a many body state (see Appendix[D] 
for details). 

We confirmed numerically that the application of 
S( g T) y ^T to the degenerate Z& state defined on the (gT) y 
torus yielded only zero energy states of the (gk + 1)- 
body model interaction. Focusing first on the densest Z& 
states, we numerically verified this property by applying 
S( g T) y ^T to the Laughlin (k = 1) and the Moore-Read 
(k = 2) states. We checked that the property was true 
at least up to TV < 13 (g = 2), N < 12 (g = 3), N < 14 
(g = 4) when km 1, and for TV < 14 when k = 2 and 
9 = 2 - 

Decreasing the filling fraction, we confirmed that the 
above property also held true for quasihole states. Again, 
we observe the completeness of the symmetrized quasi¬ 
hole subspace. We numerically verified this statement 
in a number of cases including all cases with k m 1, 
g = 2, 3, 4, TV < 6 and up to 3 added flux quanta (com¬ 
pared to the number of flux quanta for the densest 7L^ g 
state) (for N — 7, we checked all cases with one added 
flux quanta). For gk = 3, 4, we numerically verified this 
property for N < 10 and the smallest fraction of flux 
added to obtain a quasihole state. 


III. NEUTRAL EXCITATION MODES ON THE 
TORUS 

All properties discussed so far regard the construction 
of zero energy states of the (k + l)-body interaction. 
While these Z& Read-Rezayi ground states are relatively 
well understood, little is known about their neutral low 
energy excitation modes on the torus beyond the Laugh¬ 
lin case k = i) 34 l 35 l 51 l 

Most of the existing literature concentrates on the 
sphere geometry, where the neutral excitations above the 
Moore-Read and Z 3 Read-Rezayi states have been stud¬ 
ied in Refs. l2Qti22ll52l and[53]and in Ref.l23l respectively. 
In contrast, we are not aware of any related study on 
the torus geometry for the Z 2 , or for any Z& states with 







FIG. 4. Dispersion of the neutral excitation mode above 
the Moore-Read ground states on the torus as a function of 
the center of mass and relative momenta of the many-body 
states in the FQH Brillouin zone. Interpolations of the dis¬ 
persions are obtained from energy eigenvalues of finite size 
systems for a) N = — 16 and b) N — = 17. In ac¬ 

cordance with the momentum sectors Tab. [I] these systems 
feature three ground states and one ground state, respectively 
(black dots). For energies larger than E ~ 1, a continuum of 
excited states appears (not shown). The dispersions are cen¬ 
tered around the inversion symmetric momenta. The slight 
anisotropy of the dispersive part of each mode near these mo¬ 
menta is a finite size effect. We clearly observe two types of 
neutral excitations depending on the presence or absence of a 
zero energy state at the center of the dispersive part: in the 
absence of a ground state, the neutral mode has a much more 
pronounced minimum. 

k> 2. 

In this section, we start by analyzing the neutral modes 
above the Moore-Read state using exact diagonalization. 
We identify several dispersive modes (see Fig. [4] for the 
neutral excitation mode above the Moore-Read state), 
and make a connection with the neutral modes on the 
sphere. Further, we test the validity of the symmetriza- 
tion procedure beyond zero energy states. We show that 
the symmetrization construction captures the physics of 
inherent excitations of the system correctly. Symmetriza¬ 
tion of excited states above the g-layer Laughlin state will 
yield good trial states for the excitations above the Z g . 
Instead of yielding exact states, the symmetrization pro¬ 
vides a variational scheme to approximate the dispersion 
and eigenstates of the neutral excitation modes. We will 
quantitatively benchmark this construction. 


A. Neutral excitation modes from exact 
diagonalization 

The neutral excitation mode above the Laughlin state, 
the magnetoroton mode, can be understood as the dis¬ 
persing bound state of a quasiparticle and quasihole. 
Trial states belonging to the neutral mode can be ob¬ 
tained by acting on the ground state with a (lowest Lan¬ 
dau level projected) density operator. This constitutes 
the so-called single-mode approximation , wh ich has been 
verified numerically both on the sphere 3 ^ 35 * 54 * and on the 


torus geometr}EI. 

On the sphere, the Moore-Read densest ground states 
exists only for an even number of particles. There is a 
neutral mode above this state, which is also called mag¬ 
netoroton mode, because it bears strong similarities to 
the magnetoroton mode above the Laughlin state. In¬ 
deed, it is well described by the single-mode approxima¬ 
tion 54 *] Por an odd number of particles, no zero energy 
ground state of the three-body contact Hamiltonian is 
found at filling fraction v = 1 on the sphere. In spite of 
the absence of zero energy state, there is a neutral low 
energy mode, dubbed neutral fermion mode, as numeri¬ 
cally shown in Refs. 152] and E3 Both the magnetoroton 
mode and the neutral fermion mode can be viewed as 
states that minimally violate the generalized exclusion 
principle, as was shown in Ref. [54!. In the magnetoroton 
mode, the violation appears as a quasielectron-quasihole 
pair, similar to the Laughlin case. In the neutral fermion 
mode, the violation consists of a single unpaired particle 
in a background of paired particles. (See Appendix |E| for 
details.) 

On the torus, we start out by characterizing the neutral 
modes above the Moore-Read state via the exact diago¬ 
nalization of the three-body contact interaction. We find 
that both types of neutral excitation modes above the 
Moore-Read state appear irrespective of the parity of N 
(see Fig. [4| and are actually different parts of one com¬ 
mon neutral excitation mode. This neutral mode can be 
accessed in a unified way within the symmetrization con¬ 
struction, as we will show in Sec. |IIIB| by starting from 
the low-lying excitations above the Laughlin state. 

Generically, in the low energy spectrum of the (fc + 1)- 
body contact interactions, we observe patterns of low en¬ 
ergy modes, with up to four distinct dispersing branches 
that merge with the continuum of states near the four 
inversion symmetric momenta of the FQH Brillouin zone 
BZ F qh 

ko = 0, k\ = 7re x , = 7 T( 3 y, k% = 7re x T tt e y . 

(28) 

Note that fci, k 2 and k% are not valid accessible momenta 
when N is odd, but rather are defined as reference points 
in the Brillouin zone. 

On the torus, the continuous rotational symmetry of 
the quantum Hall problem is only broken by the bound¬ 
ary conditions. Still, the physics governing the neutral 
excitation modes should be dominated by shorter length 
scales comparable to the magnetic length. We can thus 
expect that the dispersion of the neutral excitation modes 
is almost rotationally symmetric and can be plotted as 
a function of a one-dimensional momentum k. To un¬ 
veil this dispersion, we have to account for the momen¬ 
tum shifts with respect to the inversion symmetric mo¬ 
menta ( |28| ) and include a geometric factor yjL x L y /N in 
order to obtain the data collapse with all system sizes 
We thus define the linearized momentum as the appro¬ 
priately rescaled minimal distance to any of the inversion 
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N = 10 ♦ 

N = 12 ♦ 

N = 14 
N = 16 ♦ 

N = 11 ♦ 

N = 13 
N = 15 + 

N = 17 


FIG. 5. Low energy spectrum of the model three-body contact interaction (Moore-Read or Z 2 Read-Rezayi state) on the 
torus for as a function of the one-dimensional momentum k defined in Eq. ( |29| ). In each graph, we show a different domain 
of the Brillouin zone, which is highlighted in grey in the inset. The black dots (respectively crosses) in the inset represent the 
Moore-Read ground s tate s for N even (respectively odd). In each domain, k is defined as the distance to the center of the 
domain (see also Eq. (29)). a) For both N even (dots) and N odd (crosses), there is a zero energy eigenstate lying at ko = 0. 
All points fall onto the same curve, defining a neutral excitation (the magnetoroton mode), b) A zero energy state is found 
at each inversion symmetric momentum in this domain for N even only. The neutral excitation energies fall onto two distinct 
curves, with the repartition determined by the parity of N. c) No zero energy ground state is found at ^3 — ( 7 r, 7 r). All neutral 
excitation energies fall onto the same curve. 


symmetric momenta 


« : = Y —^ x min { \k - ki\ | i = 0, • • • , 3} . (29) 

We show the low energy spectrum of the three-body 
contact interaction in Fig. [5] for systems with up to 
N = 17 bosons. For clarity, we represent the three in¬ 
equivalent regions of the Brillouin zone centered around 
feo, &i / &2, an d &3 in separately panels a), b) and c), 
respectively. We can clearly identify two modes with dif¬ 
ferent dispersion relations. Above the four ground states, 
there is a mode flattening out at large momentum with¬ 
out forming an energy minimum [seen in Fig. [5] a) for all 
N and in Fig. [ 5 ] b) for N even]. This mode is similar 
to the magn etoroton mode observed above the Laugh- 
lin state** 4 “ Secondly, around inversion symmetric mo¬ 
menta that do not harbor a ground state for a given N , 
there is another type of neutral mode that features a clear 
but soft energy minimum before flattening out [seen in 
Fig. [ 5 ] c) for all N and in Fig. [ 5 ] b) for N odd]. 

These characteristics of a mode without minimum and 
a mode with minimum are similar to those of the mag¬ 
netoroton and the neutral fermion modes as observed for 
the three-body contact interaction on the sphere geome¬ 
try, for an even and odd number of particles, respectively. 
Further, we observe that the first mode shows a notice¬ 
able finite size effect, while the second one is far better 
defined. This is consistent with the interpretation given 
to either mode in terms of weakly interacting elementary 
excitations on the sphere. In a background of paired 
quasiparticles, the magnetoroton mode is interpreted as 
the dispersion relation of an interacting quasiparticle- 
quasihole pair (i.e., two a quasiparticles of the under¬ 
lying Ising field theory), while the neutral fermion mode 
corresponds to the energy of an unpaired quasiparticle 
(i.e., one Ising quasiparticle). Since two quasiparticles 
induce more finite size effect in a system than one, the 


magnetoroton mode shows more finite size effects than 
the neutral fermion mode, as was observed for example 
in Ref. |52] This interpretation is relatively natural on 
the sphere, where the magnetoroton mode (respectively 
the neutral fermion mode) only appears at an even (re¬ 
spectively odd) number of particles. It is however less 
obvious that this holds on the torus geometry where both 
the magnetoroton and the neutral fermion mode are ob¬ 
served at the same filling. 


B. Approximating the neutral modes using 
symmetrization 

In Refs. [20l 1221 trial wave functions for the neutral 
modes above the Moore-Read state were obtained by 
symmetrizing over excitations of the Laughlin state on 
the sphere. If one layer is in a Laughlin ground state, 
while the other is in a magnetoroton state, one ob¬ 
tains an approximation to the Moore-Read magnetoro¬ 
ton stat d 20 * 2 ^ . On the other hand, if one layer is in a 
Laughlin quasihole state and the other is in a Laughlin 
quasielectron state, the symmetrized state approximates 
a neutral fermion state 21 . A thin torus picture provides 
an intuitive understanding of this approach, as discussed 
in App. [Ej 

To test this bilayer construction on the torus geometry, 
we have performed an extensive numerical study. We 
found that either periodic boundary conditions or twisted 
boundary conditions can be used to obtain a complete set 
of trial states for the entire neutral excitation branch of 
the Moore-Read state on the torus. We obtained good 
quantitative agreement both regarding the dispersion of 
the neutral excitation modes as well as the wave function 
overlap. 

We now compare the symmetrization construction of 
the neutral modes with the results from exact diagonal- 
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ization. 


1. Neutral excitations from the bilayer torus with periodic 
conditions 


In the following, we will see that the alternative sym- 
metrization with twisted boundary conditions developed 
in Sec. HE 2 also yields a good approximation of the neu¬ 
tral modes above the Mo ore-Re ad state. 


We start with the bosonic Moore-Read state on the 
torus with an even number of particles N. As a trial 
wave function for its neutral excitations we use 


i*r> = <s 2 - 1 (i* 1 >®i*r», (30) 


i.e., the symmetrized product of one layer with the 
Laughlin ground state |Ti) and one layer with its neutral 
excitation |Tf x ), the magnetoroton mode. 

In Refs. l20l and [2R |Ti) and |Tf x ) were obtained us¬ 
ing the composite fermion construction 4 ^. Here, we use 
the states resulting from exact diagonalization of the 
model two-body contact interaction, since the compos¬ 
ite fermion construction on the torus is at best a tedious 
task 4 ^. Note that the symmetrization procedure does not 
preserve the k x quantum number. Specific linear combi¬ 
nations of the resulting states have to be formed to ob¬ 
tain eigenstates of the corresponding translation operator 
Tf el . The overlap of the states constructed this way with 
the Moore-Read neutral excitation mode is of the order 
of 0.99 for well-defined magnetoroton states, i.e., states 
below the continuum of excitations that are found above 
a critical n. This is the same order of magnitude as the 
overlaps obtained on the sphere in Refs. [20] and HU 

When the number of particles N is odd, we use the 
following trial state 


i*r 


) — <? 2^1 ( 




qh 




(31) 


i.e., the symmetrized product of one layer with one 
Laughlin quasihole state T 


and one layer with a 


Laughlin quasielectron state |T 4e ) (obtained using ex¬ 
act diagonalization). Equation ( |3lj ) can be used to gen¬ 
erate trial states for the Moore-Read ground state |\k 2 ) 
as well as its neutral excitations |T 2 X ). Although this 
method does not yield any exact zero energy state of 
the three-body contact interaction, the approximation it 
provides for the Moore-Read ground state is fairly good 
(with an overlap of 0.998 with the Moore-Read ground 
state for N = 17). We expect this trial wave function 
to become the exact Moore-Read ground state should we 
use the exact composite fermion wave functions for T 4h 
and T^ e . This conjecture is motivated by results on the 
sphere that will be discussed in Sec. IV Likewise, T 2 X 
is a good approximation of the neutral excitation states 
(with overlaps of the order of 0.99 for states well into the 
neutral mode). The energies of the trial states for both 
N even and odd are represented in Fig. [6] along with the 
exact energies. 

Hence, this method yields very satisfactory trial states 
for the neutral excitations above the Moore-Read state 
on the torus. 


2. Neutral excitations from the bilayer torus with twisted 
boundary conditions 


As a trial wave function for the neutral excitation 
states above the Moore-Read state, we now use the fol¬ 
lowing wave function 


|^2*e) - S(2T) x ^T |^D • 


(32) 


This is the result of the action of the symmetrization 
operator on a Laughlin magnetoroton state defined on a 
twice enlarged torus (2 T) x . Similarly, applying the sym¬ 
metrization operator S( 2 t) v ^t on a Laughlin magnetoro¬ 
ton state defined on the (2 T) y torus yields another trial 
state for the Moore-Read neutral mode. 

These states have high overlap with the exact neutral 
mode states (of the order of 0.99 for N = 13). For a given 
system size, the overlaps with the trial states provided by 
the different symmetrization constructions [see Eqs. (30), 
(31), and (32)] have very similar values. We plot the en¬ 
ergies of |^ 2 ^r) i n Fig. pi (|T 2 x y ) yields the same spec¬ 
trum up to one percenr accuracy), and compare them 
to the exact spectrum of the three-body Hamiltonian. 
We note that the two subspaces created by acting on 
the Laughlin magneto-roton states with either S( 2 t) x ^t 
or are distinct, and individually not invariant 

upon a 7r/2 rotation (as already noticeable for the ground 
state). However, we numerically verified that these two 
subspaces are related by a rotation of 7r/2. 

We have thus verified that the symmetrization meth¬ 
ods provided satisfactory approximations to the neutral 
modes above the Moore-Read state on the torus. As ex¬ 
pected, the symmetrization scheme previously known on 
the sphere is equally valid on the torus, but can also 
be extended to a torus with twisted boundary condi¬ 
tions. It is however harder to reach large system sizes 
with this new scheme, because it would imply calculat¬ 
ing the Laughlin magnetoroton states for larger number 
of particles, a great numerical challenge. The same lim¬ 
itation arises when approximating the neutral mode for 
larger k : we expect the symmetrized Laughlin magne¬ 
toroton states on a g times enlarged torus to yield good 
trial states for the 7L g neutral mode. However, it is nu¬ 
merically very difficult to compute the Laughlin neutral 
mode for more than N = 13 particles. For this number 
of particles, the neutral modes above the 7L g states with 
g > 2 still show large finite size effects. 


IV. SYMMETRIZATION ON THE SPHERE 

The main result from the previous section are two 
alternative symmetrization schemes to construct exact 
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FIG. 6. Low energy spectrum of the FQH on the toms as a function of the modulus of the momentum ft, for the model 
three-body contact interaction. We compare the energy of the trial states obtained by symmetrization (hollow squares) to the 
result of exact diagonalization (filled squares and crosses). The variational states are obtained by symmetrizing a decoupled 
bilayer system (a, b, c) or a single layer system defined on the (2 T) x toms (c, d, e). The energies of these last states differ by 
less than a percent from the energies of the states obtained by the symmetrization of a system defined on the ( gT) y system. 
We thus only show the former. See Sec. EH for details on the construction of the trials states. In each graph, we show a 
different domain of the Brillouin zone, which is highlighted in grey in the inset. The black dots (respectively crosses) in the 
inset represent the Moore-Read ground states for N even (respectively odd). In each domain, ft is defined as the distance to 
the center of the domain (see also Eq. (29)). 


ground states and trial excited states of the Z g k Read- 
Rezayi type from Z& parent states on the torus geome¬ 
try. These symmetrization schemes were motivated by 
the geometrical equivalence between the double-layered 
torus with twisted boundary conditions and a twice as 
large torus. 


to a smaller single-layer sphere in this section. 


A. Z k Read-Reazayi states from symmetrization on 
the sphere 


Fractional quantum Hall states have been studied on a 
variety of other manifolds, aside from the torus. For ex¬ 
ample, model wave functions are most commonly written 
in a planar geometry, while finite-size numerical studies 
are often also carried out on the sphere. In particular, 
the multilayer construction of Read-Rezayi states was 
first numerically verified 19 ) on the sphere. While there is 
no obstruction to use the multilayer construction for any 
of the Read-Rezayi states on the sphere or on the plane, 
as there was on the torus, one can ask whether these 
geometries also allow for an alternative symmetrization 
scheme. In this section, we will show that this is indeed 
the case, i.e., the symmetrization scheme developed in 
this article for the torus can be adapted to the sphere ge¬ 
ometry or to any genus zero manifold (including the disk 
and the cylinder). Moreover, this will provide us with an¬ 
other elegant way to generate trial states for the neutral 
excitation modes in these geometries. In particular, we 
will give an alternative symmetrization construction on 
the sphere that consists in mapping a single-layer sphere 


The construction of Sec. IIB is motivated by the ge¬ 
ometrical equivalence between a ^-layered torus with 
twisted boundary conditions and a g times longer torus. 
The symmetrization is a projection on single-particle 
eigenstates with a fixed eigenvalue e 27Tls 1 9 under a trans¬ 
lation of a g -th fraction of the length of the long torus. 
While the sphere or plane does not admit an analogous 
geometrical manipulation, we can nevertheless define the 
equivalent operation to the translation on the torus. It is 
the rotation by 27 r/g, that respects the symmetry of the 
cylinder, the sphere, and the plane (disk). We are lead to 
define a symmetrization operation on rotationally sym¬ 
metric manifolds as the projection on the single-particle 
eigenstates with eigenvalue e 27Tls 1 9 under the g-th frac¬ 
tion of a full rotation (see Fig. [7]). We anticipate that this 
projection has a particularly simple representation in the 
basis of the single particle orbitals, if the chosen gauge 
is invariant under the same rotation. On such rotation- 
ally symmetric geometry, the single-particle orbitals are 
eigenstates of the angular momentum L z . Implementing 
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FIG. 7. Geometrical interpretation of the alternative sym- 
metrization prescription on rotationally symmetric manifolds, 
such as a) the cylinder or b) the disk. In either case, the 
double layer is joined into a singly connected surface via a 
topological twist defect. 


the symmetrization operation amounts to projecting on 
single-particle orbitals whose angular momentum j is a 
multiple of g. On the polynomial part of the wave func¬ 
tion, the symmetrization can be written in first quanti¬ 
zation as 


Sgs^s^ (zi, zn) = 


E * 

{0<(5i<p} 


/ 2 7T i <5 x 

[ e ^ 


1 





The sum over all 0 < Si < g with Si G Z ensures 
that only integer powers of the coordinates Zi survive in 
Sgs^s'k (zi, zn), so that the wave function is single¬ 
valued and thus physical. Using the stereographic pro¬ 
jection, one can translate this procedure from the disk 
to the sphere. The second-quantized representation of 
the symmetrization operation on the sphere or disk is 
exaclty the same as Eq. (25) on the torus. From the 
same argument that the one developed for Eq. (15) on 
the torus, we immediately deduce that Eq. (33) vanishes 
when we put (g + 1) particles at the same position if T 
is a Laughlin ground state or quasihole state. Moreover, 
the vanishing power of the symmetrized state is identi¬ 
cal to the one of T. Notice that even though we focus 
on the Zfc Read-Rezayi series in this article, our conjec¬ 
ture holds true for any (fc, r) clustered stated^ 37 . More 
specifically, Eq. (33) maps any (fc,r) clustered state onto 
a (gk, r ) clustered state. 


We now focus on the sphere geometry. There, the sym¬ 
metrization operation S g s^s precisely maps the Hilbert 
space of a sphere with gS orbitals in the lowest Lan¬ 
dau level to that of a sphere with S orbitals, once the 
correct normalization factors are put in place. A com¬ 
plication that was not present on the torus is the shift 
of the sphere, i.e., the fact that the number of orbitals 
S in the lowest Landau level is not equal to the num¬ 
ber of flux quanta N^. We start from a system with a 
number of orbitals gS = g(N ( f ) + 1) commensurate with 
g on the sphere. Acting on a many-body state with the 
symmetrization operator S g s^s divides the number of 
one-body orbitals by g. The symmetrization thus maps 


a sphere with N™ flux quanta into a sphere with 


N+ 


[ 


Nf + l 

9 


- 1 


(34) 


flux quanta, where [J is the floor function States on the 
sphere are characterized by a total angular momentum 
quantum number L z along the z-direction. For a system 
with N™ flux quanta, an eigenstate |T (L z )) is written 
as a superposition of Fock states |A) such that 


I *(£*)>= E ma>, ( 35 ) 

A) ^^7v in (Lz) 

4> 


where b\ G M and H N in (L z ) is the Hilbert space re- 

<t> 

stricted to the constraint 


iV 0 

E 0 - N72) n j (A) = L z . (36) 

3=0 

As shown in Ref. 36, the Z& states are Jack polynomi¬ 
al^^. Up to geometrical and occupation factors, their 
components in the occupation number basis are thus ra¬ 
tional numbers. We can write 


N T \f nji ] 

u tt j ’ N * 

5A = r Tl ^ttt ’ 
7=7) V n i (A)- 


(37) 


where A/^m i s the normalization factor of the orbital j, 
and r\ is a rational number. 

Due to the curvature of the sphere, single particle or¬ 
bitals on this geometry - unlike single particle orbitals on 
the torus - do not all have the same normalization factor 
for a given AU, i.e., in general, AE N in ^ A/E N in if j j- j'. 

j ’ (p j ’ <p 

To go from the large sphere to the small sphere using 
the symmetrization operator S g s^s-> one thus needs to 
adjust the normalization factors accordingly. The sym¬ 
metrized state reads 

I* (£,))= E b °( A) F(A)>, (38) 

I A) 7y in (Lz) 

4> 

where \s (A)) is defined by the occupation numbers 

n j I s (A)] = n gj+s (A), j = 0, • • • , (39) 


and the coefficients b M of the symmetrized states are just 
renormalized by their partitions’ respective amplitudes 
in the two states. 


\f n 93+s (A) 

jy, gj+a,N 4t 

1 n gj+s (^)- 

TT 9_1 \r n ai+iW \, 

L Ui =° ^gj+i,N£ V 

nf=0 n 93+i (A)! 

(40) 


Aq, 

k = E 5a II 

A:s(A)=/x j= 0 


where 0 < s < g is the sphere equivalent of the integer 
defined in Eq. (25) on the toms geometry. In terms of 
the rational coefficients, this reads 


r M = 


N+ 

E r*n 


A: s(\)=/i j =0 


n gj+s (A)- 

n?=0 n 9j+i M* 


(41) 
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Starting from a state with fixed L z , the symmetrized 
states generically have weight in different L z sectors. 


B. Numerical results 

We numerically checked that the application of the 
above described symmetrization operator S g s^s to a Z& 
state yielded zero energy states of the (gk + l)-body con¬ 
tact interaction. It is convenient to work in the rescaled 
basis where all state components are rational numbers r\. 
This allows for the exact comparison of states, despite 
the fact that we are using numerical tools. (In contrast, 
for the torus discussed in the previous section, to word¬ 
ing “exactly equal” meant always “equal up to numerical 
precision”.) We are thus effectively probing the proper¬ 
ties of some Jack polynomials under symmetrization, and 
our conclusions will apply to the Jacks themselves. 

As we are interested in finite-size systems, we have to 
account for an altered relation between the number of 
particles N needed to obtain the Z& Read-Rezayi ground 
state on a sphere with flux quanta 

kiN^ -2) = 2N (42) 

as compared to the relation on the torus Eq. ©• Given 
this constraint, the number of single particle orbitals on 
the large sphere is rarely commensurate with g. The 
polynomial nature of the Z& wave function allows us to 
easily solve this problem: adding an additional empty 
orbital at the end of an occupation number configuration 
is tantamount to multiplying each monomial by a factor 
of 1, which leaves the wave function unaltered. 

We numerically confirmed that the application of 
Sgs^s to any zero energy state of the (k + l)-body con¬ 
tact interaction yields a zero energy state of the (gk -b 1)- 
body contact interaction. Focusing on the most com¬ 
pressible state with a given number of particles, this 
property was verified for k = 1 and g = 2, 3, 4 with 
up to 12 particles, and for k = 2 and g — 2 with up 
to 12 particles. Moreover, the quasihole space built by 
symmetrizing all Z& quasihole states is complete (i.e. it 
spans the full Z^ quasihole space for the same number 
of particles and a number of single particle orbitals re¬ 
duced by one unit and divided by g). This property was 
verified for k = 1 and g = 2, 3, 4 with up to 7 particles 
and 3 added flux quanta, and for k = 2 and g = 2 with 
up to 7 particles and 3 added flux quanta. Note that in 
general, the symmetrized states are not Jacks, but some 
linear combination of Jacks. 

Finally, we used this procedure to produce trial wave 
functions for the neutral mode above the Moore-Read 
state on the sphere. The initial states are the neutral 
low energy states above the Laughlin state. These states 
can be approximated using the composite fermion con¬ 
struction-^, by considering a Laughlin system with one 
quasielectron and one quasihole. Since composite fermion 
states have exact nontrivial vanishing properties,^ we 


expect that the symmetrized state will also possess sim¬ 
ilar features. Symmetrizing the magnetoroton states, 
with g = 2, yields one state per value of L, except for 
L — 0 and L — 1 if N is even (respectively L = 1/2 if N 
is odd). Note that there is no neutral mode state in these 
sectors. For the accessible system sizes (up to TV = 14 
particles), we also verified that the trial states, although 
different from the ones obtained in Refs. I2QH221 provide 
an equally good approximation to the neutral modes. 
The symmetrization construction described in this sec¬ 
tion thus provides a unique method to approximate both 
the magnetoroton and the neutral fermion modes on the 
sphere. In contrast, the methods of Refs. [2Qlf22l all used 
two different sets of initial states to obtain these modes. 

We have thus provided an alternative symmetrization 
method on the sphere. Z^ states can be obtained by 
symmetrizing a Z& state defined on a g times larger 
sphere. For a state written in the occupation number 
basis, the symmetrization procedure consists in select¬ 
ing one in g orbitals and discarding all of the other 
ones. This procedure was directly adapted from a simi¬ 
lar symmetrization scheme on the torus, accounting for 
the sphere specificity in the normalization of the single 
particle orbitals. We numerically confirmed that all Z^ 
quasihole states could be obtained using this method. 
This highlights a previously unknown property of some 
of the Jack polynomials. In contrast, the alternative sym¬ 
metrization procedure that could be derived from the 
second quantized representation of S^ q t) x ^t ( see Ap¬ 
pendix [C]) destroys the squeezing hierarchy of the occu¬ 
pation number basis on the sphere, and therefore cannot 
be used. 


V. CONCLUSION 

This study was motivated by the fact that the previ¬ 
ously known multi-layer symmetrization construction of 
Z fe Read-Rezayi states misses several zero energy state 
of the model Hamiltonian on the torus. Our result is 
an alternative projective construction scheme that turns 
out to be as powerful as the multi-layer symmetrization 
construction and ideally complements it. This novel con¬ 
struction obtains a 7L g k state from a Z& by reducing the 
size of the manifold - torus or sphere - on which the 
state is defined by a factor g. The construction has a 
suggestive geometrical interpretation in terms of topo¬ 
logical twist defects connecting the different layers. On 
the sphere or the plane, it leads to a conjecture regarding 
the mathematical properties of some Jack polynomials, 
and more generically the clustered states. 

Beyond exact statements on zero-energy states, we 
showed that the projective constructions also yield excel¬ 
lent approximations to the collective low energy neutral 
excitation modes above the Z& Read-Rezayi states. 

Our results open several natural directions for future 
work, such as the extension to other fractional quantum 
Hall states, in particular those of fermions. 


14 


ACKNOWLEDGEMENT 


We acknowledge Z. Papic, B. Estienne and M. 
Barkeshli for discussions. We are especially grateful 
to Curt von Keyserlingk for helpful comments. The 
authors acknowledge financial support from DARPA 
SPAWARSYSCEN Pacific N66001-11-1-4110. B.A.B 
and N.R. were supported by NSF CAREER DMR- 
095242, ONR-N00014-11-1-0635,MURI-130- 6082, MER- 
SEC Grant, Packard Foundation, Keck grant and the 
Princeton Global Scholarship. NR and CR were sup¬ 
ported by ANR-12-BS04-0002-02. 


Appendix A: Translation operators and momentum 
quantum numbers 


We consider a rectangular torus spanned by the vectors 
L x e x and L y e y , where e x and e y are two unit vectors. 
Translation operators on the torus can be factorized into 
the product of a center of mass and a relative transla¬ 
tion. The center of mass translation operator along the y 
axis and the relative translation operator along the x axis 
commute with each other and with the Hamiltonian. The 
eigenstates of the Hamiltonian thus carry the correspond¬ 
ing momentum quantum numbers k that belong to the 

FQH Brillouin zone BZfqh = jfe = ^k x e x + |^k j 

with k x = 0, • • • , GCD(7V, N^) — 1 and k y = 0, • • • , N^ — l. 
Here, GCD stands for the greatest common divisor, and 
e x , e y are such that e* • e.j = Sij. Focusing on the densest 
Zfc states, for k odd (respectively even), the momentum 
quantum numbers k x ,k y thus belong to a N^/2 x 
(respectively N ^ x N^) Brillouin zone. These momen¬ 
tum quantum numbers are defined in the Landau x-gauge 
with vector potential A(r) = (0, —Bx). 

Due to the center of mass symmetry, there is a q- 
fold degeneracy of all states in the FQH Brillouin zone, 
where q = A^/GCD(iV, N^). Without losing informa¬ 
tion, we can thus work in a reduced Brillouin zone 57 of 
size GCD(7V, N^) x GCD (A, 7V 0 ) 


-oyred 

^FQH 


| . 27r 27T 

i k — — k x e x H - — k v e 


y^y 


k* = <),••• ,GCD(7V, Nfy) — 1; 
ky = 0, • • • ,GCD(7V,7V 0 )-lJ. 


(Al) 


In the case of the Z& Read-Rezayi densest state, the pre¬ 
vious equation leads to a reduced Brillouin zone BZpg H 
of size N^/2 x N^/2 (respectively x N^) for k odd 
(respectively k even). Note that for k even, the reduced 
Brillouin zone coincides with the full Brillouin zone, while 
for k odd, BZpg H corresponds to half of the full Brillouin 
zone. 



FIG. 8. Sketch of the two possible symmetrizations for the 
double-layer toms with twisted boundary conditions in both 
directions. We look at the simple case where L x = L y . (a) 
and (d) show the initial system where the boundaries have 
to be glued according to the plain, dotted or dashed lines at 
the edge, (a) requires the equivalent single-layer system (b) 
to be defined by the spanning vectors L x e x + L y e y of length 
y/L% + L| and 2 L y e y . Analogously, the single-layer system 
(e) for (d) is set by the spanning vectors 2 L x e x and L x e x + 
L y e y of length + L|. From (b) [resp. (e)], we apply 

the symmetrization procedure S{ 2 t) v ^t (resp. S( 2 t) x ->t) to 
obtain a Moore-Read state (c) [resp. (f)]. In both cases, one 
ends up with a Moore-Read state on a twisted torus. Notice 
that the orientation of the orbitals (along L 2 ) depicted here 
is chosen accordingly to the gauge used when expending the 
wave functions on the many-body basis. 


Appendix B: Symmetrization on a double-layer torus 
with twisted boundary conditions in both directions 


In this Appendix, we discuss the result of a sym¬ 
metrization over a multi-layered torus with twisted 
boundary conditions in both the x and the y directions. 
This manifold is geometrically equivalent to a single-layer 
torus as well. For simplicity, let us consider a double¬ 
layered torus spanned by two orthogonal vectors L x e x 
and L y e y . It can be related to a single-layer torus in two 
ways (see Fig. |8|. 


One option is to account for the twisted boundary con¬ 
ditions in x-direction by taking one spanning vector of 
the single layer torus as Li = 2 L x e x . Then, to also ac¬ 
count for the twisted boundary conditions in the y di¬ 
rection, the other spanning vector must be chosen as 
L 2 = L x e x + L y e y . Thus, a rectangular double layer 
torus spanned by L x e x and L y e y with twisted bound¬ 
ary conditions in both directions is equivalent to a single 
layer torus spanned by L\ and . By the arguments 
in Sec. II C[ this geometrical equivalence also holds on 
the level on the Hamiltonian. Symmetrization reduces 
this (£ 1 ,^ 2 ) torus to a torus spanned by L\j2 and L 2 . 
In the x-Landau gauge, the action of the symmetrization 
operator is given by Eq. (24) in this case, despite the fact 
that the torus is nonrectangular. We numerically checked 
that this symmetrization of the v = 1/2 Laughlin state 
on the (^ 1 ,^ 2 ) torus yielded the Moore-Read state on 
the (£q/2, L 2 ) torus for systems with up to 12 particles. 
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Alternatively, to first account for the twisted boundary 
conditions in the y direction, one may take one spanning 
vector of the single layer torus to be L 2 = 2 L y e y . Then, 
to also account for the twisted boundary conditions in x 
direction, the other spanning vector must be chosen as 
L\ = L x e x + L y e y . Again, a rectangular double layer 
torus spanned by L x e x and L y e y with twisted boundary 
conditions in both directions is equivalent to a single layer 
torus spanned by these vectors L\ and L 2 . By the argu¬ 
ments in Sec. |II C[ this geometrical equivalence also holds 
on the level on the Hamiltonian. Symmetrization reduces 
this (Li,L 2 ) toms to a torus spanned by L\ and T 2 /2. 
In the x-Landau gauge, the action of the symmetrization 
operator is given by Eq. (25) in this case, despite the fact 
that the toms is nonrectangular. We numerically checked 
that this symmetrization of the v = 1/2 Laughlin state 
on the (Li,L 2 ) toms yielded the Moore-Read state on 
the (Li, L 2 / 2 ) toms for systems with up to 12 particles. 


Appendix C: Practical implementation of S^ g T) i ^T 

In this Appendix, we provide some technical informa¬ 
tion to implement both S( 9 t) x ^t and S( g T) y ^T- As we 
will show, their action is simple to write once the wave 
function is decomposed on the occupation basis. 

Enlarging the torus length in the x-direction — We 
detail this procedure sketched in Fig. [ 3 ] c) and Eq. (24) 
in the case, where we symmetrize a Z& state on a toms 
of size ( gT) x , pierced by gN^ flux quanta to obtain a 
candidate for a 7L g & state on a T toms pierced by flux 
quanta, by projecting on orbitals with quantum numbers 
j that are integer multiples of g. For a system with gN^ 
flux quanta, a state |T (k^)) in a given k y sector is written 
as a superposition of Fock states |A) such that 


l*(M>= E |A), 


(Cl) 


Ityt'UgN'f, (k y ) 


where b\ G C and HgN# (k y ) is the Hilbert space re¬ 
stricted to the total momentum constraint 


( gN^ — l 


K = E mod > 


(C2) 


3=0 


with k y G {0, • • • , gN$ — 1}. Here and below, a Fock state 
| A) corresponds to the occupation-number configuration 
n (A) with 


n (A) = {nj (A), j = 0, • • • , — 1} , (C3) 

where nj (A) G No is the occupation number of the single¬ 
particle orbital with momentum j. 

The symmetrized state S^ 9 t) x ^t (k y )) appears in 
the k y mod momentum sector and reads 

S (gT) ^ T |¥ (k y )) = E I s * ’ ( C4 ) 

|A)G'HyiV </) (ky) 


where | s x (A)) is defined by the occupation numbers 
9 - 1 

n j I s * (A)] = E n j+iNt (A), j = 0, • • • , Nj, - 1 (C5) 

i =0 

and the coefficients b y of the symmetrized states are given 
by 


n a -i 


b [ZU rij+iN,, (A)_ 

^ bx n niA^(A)! 


A : s x (\)=fi j =0 


(C6) 


Enlarging the torus length in the y-direction — We now 
detail the procedure sketched in Fig. [ 3 ] b) and Eq. (25) 
in the case where we symmetrize a Z& state on a toms 
of size L x x gL y pierced by gN^ flux quanta. For each 
group of g consecutive orbitals, the symmetrization op¬ 
erator selects the 8 th orbital of them (where s < g is a 
nonnegative integer), and discards all other g — 1 orbitals. 
The symmetrized state S^ g T) y ^r (k y )) is given by 

‘V)^t|*(M)= E \,(a)MA)>, (C7) 

|a )e'Ug N(l) 

where \s y (A)) is defined by the occupation numbers 

n j ( s y (^)) = n gj+s (^) 5 3 = 0 ? ’ ’ ’ 5 ^<f> ~ f (C 8 ) 

and the coefficients b g of the symmetrized states are 
N+-1 


k = e b * n 

A : s y (X)=/i j =0 


n gj+s (A)j 

n^o 1 n gj+i (A)! 


(C9) 


As a consequence of the second part of Eq. (25), all state 
components containing at least one j 7 ^ gj + s such that 
n~- (A) 7 ^ 0, are discarded. The number of particles N is 
thus conserved by S(gT) y T. 


Appendix D: Symmetrization and many-body 
momentum quantum numbers 


We denote by T x el the many-body operator that cre¬ 
ates a relative translation by L x e x /N in position space. 
The properties of the momentum quantum numbers k x 
and k y under the symmetrization S^ g T) y ^r are relatively 
simple. The former momentum is conserved, and the 
latter is mapped from k y to 

k / = K~ sN 

y 9 

where s is the integer defined in Eq. ( [25] ) (0 < s < g). 

The properties of the momentum quantum numbers 
under the action of S^ q t) x ^t are slightly more compli¬ 
cated. We detail these properties in the following para¬ 
graph. 
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When the symmetrization operation defined in 
Eq. (25) is applied to a state on the (gT) x toms with cen¬ 
ter of mass momentum k y mod (gTV^), it yields a state 
on the T torus with the momentum k y mod N^. The 
relative momentum quantum number k x is not generally 
conserved during this transformtion. However, in some 
cases, the symmetrized state is still an eigenvector of the 
relative translation operator. We examplify this feature 
on the case of k = 1 (i.e., the initial state is the Laugh- 
lin state), for a symmetrized state with filling fraction 
v = g/2 , and compute the corresponding eigenvalue. We 


will treat the cases of even and odd g separately. 

Case g odd — We start from a state on the (gT) x 
torus |T (k y )). Since the Brillouin zones of both the orig¬ 
inal and the symmetrized tori are twice as large as their 
reduced Brillouin zone, a state on either of these systems 
can be an eigenstate of (T x el ) 9 . On the (gT) x torus, an 
eigenstate of (T x el ) 9 with relative momentum k x writes 

^-i 

27ripk cc 

\4>(K,ky))= £ eW (7V e1 )^ l^(k y ))- (D2) 

p =0 


Note that at filling fraction g/2 with g odd, must be even. The action of the symmetrization operator on this 
momentum eigenstate is given by 




-1 


S {gT) x =>T M) “ £ e aN ’t> /2 (rr') 9P S(gT) x ^T l’J'(M) 

p=0 

g-lN^/2-1 2 ,i(j,+jAr^/2)k a , , , . ;o> 

= ^2 £ e-- (l-rel\9(P+jK4>/2) 

3=0 p =0 


(7r') £ 


S(gT) x =>T I ^)) 


9~ 1 

E 

3=0 


27rijk x 

e a 


27rip(k x /g) 


NJ 2-1 

£ e^^(r x °y p s {gT) ^ T Mk y )), 

p =0 


(D3) 


where we have used the fact that (7 ~^ ei y N(f> ^ 2 = (T x el ) N = 1. If k x is not a multiple of g , the result is zero due to the 
factor Yl 9 j =o e 9 > and there is no symmetrized state. If k x is a multiple of g , however, the symmetrized state is an 

eigenstate of (T x el ) 9 and the corresponding eigenvalue is 


x g' 


(D4) 


Case g even — The filling fraction v = g/2 of the symmetrized state is an integer, therefore can have either 
parity. When is even (i.e., TV is a multiple of g), the symmetrized state is not an eigenstate of the relative 
translation operator. We focus on the case where is odd. Due to their respective filling fractions 1/2 and g/2 E N, 

the original and symmetrized states may be eigenstates of ( T x el ) 9 and (T x el ) 9 ^ 2 , respectively. 


We start from a state on the (gT) x torus, with relative momentum k x , such as defined in Eq. (D2). The symmetrized 
state writes 


s£> _ 


271 -ipka; 


S(gT) x ^T \4>{k x , ky)) — £ e N * ( - 9/2) (JZ* A ) 9P S(gT) x ^T |^(ky)) 

p =0 

9/ 2—1 N ^~ 1 2ni(p+jN (f) )k x ( , • -\ T \ 

= E £e ^ (9/2) * S(gT)^ T \*{k y )) 


3=0 P=0 

S/2-1 

E 27r ij’ k a 

e «/=. 

3=0 


N *-1 


2-K\p(2k x /g~) 


£ e ^ (^y p s igT) ^ T mk y )), 

p =o 


(D5) 


where we have used the fact that (' T x el ) dN(p = (' T x el ) 2N = 1. If k^ is not a multiple of g/2 , the symmetrized state 
vanishes due to the Fourier sum. Otherwise, the above expression can be reexpressed in the following way (up to an 
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overall factor gj 2) 


2 7 rip(k cc /g) ^ 1 27rip(k cc /g) 

S(gT) x ^T\<t>(kx, kj/)) — X] e W<#,/2 (7T') e ^ /2 (7^ el ) 9P 5( ffT ) x _,. T ^(k^)) 

p=0 P=(N<t,+ 1)/2 

JV*-1 

J—^ 27riq(k a; /g) a 

= E e ^ (TTr^T^Tl^M) 

q=0, even 

Ar ^'“ 1 2,ri(q + iV </) )(k a ./ fl ) g g 

+ E e (7T 1 ) ‘«V ) .-r|*(k l ,)> 


q= 1, odd 

iV ^ _1 27riq(iV 0 + l)(k x /g) g 

= E e dr 1 ) 59 

< 2=0 


S(gT) x ^>T WM) ’ 


(D6) 


The symmetrized state is thus an eigenstate of (TJ el ) s ^ 2 with eigenvalue 

V! x = + 1 ) kg m od . 


(D7) 


Appendix E: Thin torus perspective on 
symmetrization 

1. Thin torus and zero energy states 

To intuitively understand which states appear under 
the symmetrization operation, we can consider the action 
of symmetrization on the root configurations or on the 
thin torus limit configurations that represent many-body 
states. Theses are single occupation number configura¬ 
tions in the single particle basis labelled by linearized 
momentum index j on the sphere and torus geometry, 
respectively. 

On the sphere, the approach of Ref. |36] makes use of 
the clustering properties of the states to rigorously derive 
fractional quantum Hall ground states^ and quasihole 
stated^ starting from a root configuration that obeys a 
generalized exclusion principle 59 . Such a derivation of 
the state from a root configuration is not possible on the 
torus. Nevertheless, for many FQH states (in particular 
for the Read-Rezayi series), the global counting of ground 
states and more generally of zero modes (i.e., quasihole 
states) is correctly predicted^ by counting single occu¬ 
pation number configurations that obey the generalized 
exclusion principle on the torus. The physical signifi¬ 
cance of the single occupation number configurations on 
the torus is that they correspond to a charge density wave 
pattern in the limit where the circumference L x is taken 
to zercPH^. They are thus referred to as “thin torus” 
configurations. 

Let us for concreteness consider the Laughlin (k = 1) 
and Moore-Read (k = 2) states at filling v = 1/2 and v = 
2/2, respectively. Their all owed root configurations obey 
a (&, 2)-exclusion principle 36 37 , meaning that no more 


than k particles may be found in 2 consecutive orbitals. 
Periodic boundary conditions on the torus imply that the 
(&, 2) exclusion principle has to be satisfied also between 
the first and the last orbital. Imposing this constraint 
on the thin torus configurations immediately delivers the 
k + 1 degenerate ground states on the torus when the 
number of bosons A is a multiple of k. For example, the 
two topologically degenerate Laughlin states {k = 1) on a 
torus with = 6 orbitals are represented by the unique 
two densest configurations that obey this principle 

101010, 010101. (El) 

(The notation for the thin torus configurations lists the 
sequence of occupation numbers rij of the consecutive or¬ 
bitals j = 0, • • • , — 1. For the first example no = 1, 

n\ — 1, etc.) The symmetrization of the root configura¬ 
tion of two layers of Laughlin states with respect to the 
layer quantum number of the orbitals, represented by the 
operator £ 2 ^ 1 , yields 


( ioioio\ 

^ 101010 J 


= 202020 , 


(E2) 


which is one of the correct thin torus configurations of 
the Moore-Read state (k = 2). The other Moore-Read 
ground states 


111111 , 


020202 


(E3) 


can be obtained analogously by symmetrizing over differ¬ 
ent combinations of the two Laughlin ground states (El) 
in the two layers in accordance with what is listed in 
Tab. [TJ That symmetrization of the entire Laughlin state 
(not only its thin torus configuration) exactly yields the 
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FIG. 9. a) Schematic representation of the impossibility to 
find a thin-torus configuration respecting the (1,2) exclusion 
principle for a bilayer torus with N = 7 bosons, and m 7 
flux quanta, b) Inserting a topological defect that connects 
the two tori relieves the frustration. A configuration respect¬ 
ing the exclusion principle is then possible. 


Moore-Read state is a highly nontrivial fact which can¬ 
not be deduced from this observation on the thin torus 
configurations alone. However, if the symmetrization al¬ 
ready fails on the level of thin torus configurations, it 
cannot work for the entire state either. 

We can use this fact to illustrate why symmetrizing 
over several layers of the Laughlin state cannot yield all 
ground states of the Moore-Read type, or more generally 
of the Zfc Read-Rezayi type, on the torus. Turning to 
the torus with odd number of flux quanta, e.g. = 7, 
a single Moore-Read ground state exists with thin torus 
configuration 


1111111 (E4) 

that obeys the (2,2) exclusion principle. Form the thin 
torus configurations, one can understand that this state 
cannot be obtained by symmetrization over two indepen¬ 
dent layers of Laughlin tori. The problem is that there 
exists no zero energy state that obeys the (1,2) exclu¬ 
sion principle on a torus with odd and the desired 
number of particles. One would need the thin torus con¬ 
figurations to be 


<$2—^1 


0101010' 

, 1010101 , 


= 1111111 . 


(E5) 


However, due to the periodic boundary conditions, the 
lower of the two configurations, 1010101 is not allowed, 
for it violates the (1,2) exclusion principle between its 
first and last orbital (see Fig. |9|. The twisted boundary 
conditions are exactly what is needed to make t he t hin 
torus configuration on the lefthand side of Eq. (E5) al¬ 


lowed: With them in place, the (1,2) exclusion principle 
has to be satisfied between the last orbital in the upper 
layer and the first orbital in the lower layer as well as 
between the first orbital in the upper layer and the last 
orbital in the lower layer. 

Equivalently, we can apply the two symmetrization 
schemes for the (gT) y torus and for the (gT) x that we 


summarized in Eqs. (24) and (|25j) , respectively, 


1111111 = %r )y ^ T (10101010101010) 

= S (gT) ^ T (10101010101010). iEG) 
Again, the fact that these equalities hold not only for 
the thin-torus configuration, but for the entire Laughlin 
states, is a highly non-trivial fact. 


2. Thin torus perspective on neutral modes 


In terms of thin torus configurations, states belonging 
to the neutral modes appear as minimal violations of the 
generalized exclusion principle at fixed particle density. 
For example, a thin-torus configuration belonging to the 
magnetoroton mode of the Laughlin state is given by 

101101010010, (E7) 

• o 


where dot (•) and circle (o) mark the positions of 
the quasiparticle (violating the exclusion principle) and 
quasihole, respectively. Using this construction, explicit 
wave functions were derived in Ref. [54| thus validating 
this interpretation of the neutral mode. 

As discussed in Sec. |IIIB| neutral excitations above the 
Moore-Read state can be captured by a bilayer construc¬ 
tion involving symmetrization over Laughlin charged or 
neutral excitations. Let us illustrate this scheme on the 
torus geometry with the help of thin torus configurations. 
For that, we consider the symmetrization over double 
layer configurations with a single violation of the (1,2) 
exclusion principle to construct a single violation of the 
(2, 2) exclusion principle. We have two choices to ob¬ 
tain the same symmetrized root configuration: In the 
first case, one Laughlin layer is in a magnetoroton state 
and the other layer is in the Laughlin ground state, for 
example 


/ioiioioiookA 

<S 2 ^i • ° = 111202020111. (E8) 

\oioioioioioiy 

In the second case, one layer is in a Laughlin quasielec¬ 
tron state and the other layer is in a Laughlin quasihole 
state 


<S 2 -^i 


/101101010101 
1 010101010010 


= 111202020111 . 

• o 


(E9) 


Observe that Eq. (E8) is an allowed (minimally vio¬ 
lating) con figu ration for periodic boundary conditions, 
while Eq. (E9) is allowed for twisted boundary condi¬ 
tions. 
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